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Abstract 

In this paper we carry out Quantum Monte Carlo simulations of a quantum 
particle in a one-dimensional random potential (plus a fixed harmonic poten- 
tial) at a finite temperature. This is the simplest model of an interface in 
a disordered medium and may also pertain to an electron in a dirty metal. 
We compare with previous analytical results, and also derive an expression 
for the sample to sample fiuctuations of the mean square displacement from 
the origin which is a measure of the glassiness of the system. This quantity 
as well as the mean square displacement of the particle are measured in the 
simulation. The similarity to the quantum spin glass in a transverse field is 

noted. The effect of quantum fluctuations on the glassy behavior is discussed. 
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I. INTRODUCTION 



The problem of a quantum spin glass in a transverse field was recently in the center 
of theoretical and experimental interest In particular the question of the interplay 

between glassy behavior and quantum fiuctuations and the properties of the quantum tran- 
sition at zero temperature were the subject of several investigations, both for the case of the 
infinite-ranged spin glass and more realistic three dimensional models. 

A simpler model which catches many of the essential features of the (classical) spin 
glass problem is that of a directed manifold in a random medium |]5|-pil|]. An even higher 
simplification occurs for the case of a zero dimensional manifold, which is equivalent to a 
particle in a random potential (which is localized by an additional fixed harmonic potential) 
||TTHl6|. This model has been found to require a (infinite-step) Parisi type solution when 
the random potential has long range correlations [0. For a random potential with short 
ranged correlations a one-step replica-symmetry-breaking (RSB) solution has been found 
P Jl5|JT6|| . A single particle in one dimension does not have a sharp transition into a glassy 
phase. But in infinite dimensions it does. It turns out that the analytical solution which 
utilizes the variational approximation, still possess a sharp transition at finite dimensions 
(including one dimension). This occurs since the replica symmetric (RS) is not able to 
capture the glassy features of the systems once they become strong enough, and gives rise 
to unphysical results, like a non-monotonic mean square displacement of the particle from 
the origin as a function of the temperature. Below a certain temperature the RSB solution 
yields a much better physical result, and in particular the correct non-analytic behavior at 
T = as a function of the strength of the random potential in agreement with an Imry-Ma 



type argument One should notice though, that the transition which has been found 

for a particle in a random potential is of the Almeida-Thouless type ||17| , p!8| in the sense 
that it is associated with RS- RSB transition but not with an order-disorder (spin-glass like) 
transition. 

The success of the variational treatment of a classical particle in a random potential, led 
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one of us recently [lT9| , pO| to investigate a quantum analog, i.e. to turn on h and consider the 



effect of quantum fluctuations, e.g. tunneling on the glassy behavior of the particle. This 
was suggested by the recent theoretical treatment of the quantum spin glass in a transverse 
field ii- 

Analytically we have found a glassy phase characterized by RSB, which is destroyed by 
the quantum fluctuations for strong enough or alternatively for small enough m, which is 
the particle's mass. The variable Yi' jm plays the role of the transverse field in the quantum 
spin glass problem. A schematic phase diagram is depicted in Fig. 1. The full details of the 
analytical investigation are given in ||2^ , to be referred in the sequel as I. 

Recently ||21[ one of us used the model of a particle in two spatial dimensions under 
the influence of an harmonic and a quenched random potential to describe the melting 
transition of the flux lattice in high-temperature superconductors with columnar disorder. 
The so called cage model has been introduced originally by Nelson and Vinokur In 



this model a single flux line is represented by the world line of a quantum particle. The 
influence of neighboring flux lines is taken effectively as the cage harmonic potential. The 
magnitude of Ti is determined by the size of the system along the z-axis and the value of 
the temperature. This shows the usefulness of the toy-model to other physical systems of 
interest. 

Our aim in this paper is two-fold. First we carried out a quantum Monte Carlo simulation 
of the system in one dimension in order to compare with the analytical results obtained in I 
for the mean square displacement. Second, since both the random potential and the quantum 
fluctuations increase the mean square displacement, this quantity by itself is not enough to 
give a clear picture concerning the strength of the glassy behavior of the system. Hence 
we have measured in the simulation, and also calculated analytically the sample to sample 
fluctuation of (x^), which shows that the glassy behavior of the system diminishes as the 
quantum fluctuation increase until the eventual demise of the glassy phase for strong enough 
fi^ /m. This decline in the sample to sample fluctuations is a gradual effect, which analytically 
culminates in the transition from an RSB solution to a RS solution. As mentioned previously. 



in a simulation which is carried out in one dimension, we do not expect to observe any sharp 
transition. The trapping of the particle in deep local minima of the random potential gives 
rise to a sticky behavior- i.e. a freezing of the mean square displacement from the origin. 
This effect is countered by tunneling among the different minima which enables the particle 
to escape from a local minima and thus diminishes the glassy behavior. 

In the next section we define the model. In Section 3 we review some of the results 
obtained in I and obtain a new analytic result for the sample to sample fluctuations of the 
mean square displacement of the particle from the origin. In Sect. 4, we describe the details 
of the quantum Monte Carlo (QMC) simulation. Section 5 is devoted to a discussion of the 
results and a summary. A comparison is made between the simulation and the theoretical 
results obtained in I and in Section 3. In the Appendix we give some further details of the 
calculation presented in Section 3. 



II. THE MODEL 

The Partition function for a particle at finite temperature T = 1/^^/5, subject to a 
harmonic potential and a random potential V, is given by the functional integral : 

rU 



Z{U) = f [dx] exp / 

A(o)=x((7) I n Jo 



m±{uf , ^^^{uY ^ ^^^^^^^ 



du), (2.1) 
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where x is a A^-dimensional vector (A^ is the number of spatial dimensions), and U = (3h. 
The variable u has dimensions of time and is often referred to as the Trotter dimension. 
We observe that the trajectory x(m) forms a closed path. In this paper we are concerned 
with a random quenched potential V{x), which is Gaussian distributed. This means that 
the probability for a given realization of the potential is given by: 

P{y{x)) = Cexp{- J dxdx'V{x)A{x - x')V{x')), (2.2) 

with some known function A{x — x'). It is thus sufficient to know only the first two moments 
of the distribution vis. 
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(\/(x))^ = 0, (\/(x)y(x'))^ = -Nf (^^4r^) • ^^-^^ 

where the functions / and A are related to each other. The function / describes the 
correlations of the random potential. In this paper we consider two cases. One, which we 
call the case of long ranged correlations of the potential, for which / is taken to decay as a 
power at large distances: 

m = ^(Y^ («o + y) (2.4) 

with 7 = 1/2 . This case corresponds in one dimension to an interface in the random field 
Ising model P, p^|T5| . The parameter oq plays the role of a short- distance regulator for 



/. Another type of random potential we consider has Gaussian correlations and we refer to 
it as the case of short-ranged correlations. For this case, the function / is taken to decay as 
a exponential function at large distances: 

f{y) = I exp(-iy) (2.5) 

In the classical case it has been shown that (for = 1) even when correlations fall expo- 
nentially fast, the physics is equivalent within the variational approximation to the case of 
random potential with power law correlations (eq. ( p.4|) ) and 7 = 3/2 at large distances 
PJl5[| . This fact also holds in the quantum has been demonstrated in I. Thus we will 



compare the results obtained in the QMC for the distribution (|2.5|) with results obtained in 
I for the case of 7 = 3/2. 



III. THEORETICAL CONSIDERATIONS 

In paper I we investigated the problem using the replica method and the variational 
approximation. Here we review briefly some of the formalism and discuss the new derivation 
of the sample to sample fluctuation of the mean square displacement that is a measure of the 
glassiness of the system and is needed to compare with the results of the QMC simulations. 
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Some details are deferred to the Appendix. Readers who are interested only in the details 
of the QMC simulation can skip this section. 

In the replica method variational approach the system is represented by an n-body 
variational Hamiltonian: 

1 



J- /" r 1 

2 Jo V ^ J 

1 f'^ [V 

~ 7^ / du'y^Sab{u-u')Xa{u) -Xbiu'). (3.1) 

Zn Jo Jo ^ 

The matrix Sab{u — u') is determined by extremizing the variational free-energy which is 
given by: 

n/3 (F)^ IN = {Hn - hn)hjn - In / [dx]e-'^"/^ (3.2) 

Here Hn is the exact n-body Hamiltonian. The hmit n — > has to be taken at the end. 
The propagator associated with hn is given in frequency space : 

Gabiuj) = ([(ma;2 + /.)!- s(a;)]-i)„,. (3.3) 

CO is the frequency conjugate to the Trotter time variable u, and takes the values: 

^i = jjl. / = 0,±1,±2,--- , (3.4) 

and the matrix Sabi'jj) is related to Sab{u) by: 

Sab(C) = ^ exp(-ia;, C) Sab{uJi). (3.5) 

1^ l=-oo 

We have found a self-consistent solution to the variational equations where only the 
diagonal elements of the matrix Sab are 'time' dependent, and the off-diagonal elements are 
independent of the Trotter time. Thus 

Saa{uj) = Sd{uj), (3.6) 

Sab{uJ, z) ^ s{z) d^fl , a^b, (3.7) 
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where the Parisi parameter < z < 1 labels the "distance" between rephcas indices ab. A 
similar behavior follows for the propagator matrix Gah{i^) with a similar notation Gcl{oj) and 
G{z) (These are the same as the quantities raiuj) and r{z) used in I) . 
The mean square displacement from the origin is given by: 

1 oo 

{{^'))r/N^- E GaM- (3-8) 

k=—oo 

This quantity was evaluated in I. 

The sample to sample fluctuation of (x^) is a measure of the glassiness of the system. 
In the replica approach together with the variational approximation, this quantity is repre- 
sented by 

((-f ), - - i-i^k - « p-^) 

where a and b are indices for replicas. Follow the notation we used in I, the sample to sample 
fluctuation of the mean square displacement becomes: 

For a particle in N dimensions it is self averaging in the large N hmit but not ior N — 1. 
The degree of non-self-averaging is a measure of the glassy behavior of the system. 

In the Appendix we give details of the numerical evaluation of the sample to sample 
fluctuation of ((x^)) both in the replica symmetric and in the RSB phases. The calculation 
is done both for the case of continuous RSB which occurs for a random potential with long 
ranged correlations, and for the case of short-ranged correlated potential where there is a 
one-step RSB. The results are depicted in Fig. 2 and Fig. 3 respectively. We observe 
that the glassiness of the system increases with decreasing temperature, but decreases with 
increasing h^/m = 1/k. Recall that the transition temperature Tc{fi) between the RS and 
RSB phases decrease with decreasing k as was obtained in I. 

Another quantity that could also be used as a measure for the glassiness of the system, 
but we did not measure in the QMC simulation is the sample to sample fluctuation of the 
susceptibility, 
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(3.11) 



which is given by: 



R 



ix) 



R 



3/32 



1 + 



N 



dzG\z) 



- (^J^ dzG{: 



(3.12) 



IV. NUMERICAL SIMULATION 



We apphed the path integral Monte Carlo method (PIMC) in one space dimension to 
calculate the relevant physical quantities we are interested in. The partition function of this 
system at a given temperature is given by 



Z{T) 



x(C7)=x' f 1 /-^ 

[rfx] exp <^ -- / 
x(o)=x I n Jo 



du 



(4.1) 



We discretize x(r) into M+1 points with x{l) equal to x{M + 1). The partition function 
now becomes an M-dimensional ordinary integral, i.e., 



z{T) = / n 



M 



dxi 



■ exp 



1=1 



M 



(4.2) 



where e — U/M, A — ^J^mhejm. In this way we are able to apply the Metropolis method 
to perform the integration for the partition function. 

In this form the partition function of a quantum particle is similar to the classical Boltz- 
mann distribution of a polymer ring with M beads under the influence of an applied harmonic 
potential and a random potential. The beads have a harmonic spring interaction between 
neighboring beads and in addition each bead fccis a combination of a harmonic (w.r.t the 
origin ) and a random local potentials. This gives us an intuitive picture of the simulation. 
In the simulation one attempts to move each bead in its own turn and one checks if the 
"energy" ( minus the argument of the exponential) decreases or increases. Then one ac- 
tually moves the bead in accordance to a detailed balance algorithm, e.g. the Metropolis 
algorithm. 
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The problem in doing PIMC comes from the fact that for small e (or large M) the beads 
are not easy to move due to a very large spring constant, thus the acceptance rate is low for 
a reasonable size move and convergence is not easy to achieve in a reasonable time. For this 



reason many efforts have been made to circumvent the problem [Q, |3y], |3T|] . We use 
the normal mode PIMC to perform the calculation ||2^. In our program the motion of the 
'beads' includes 2 parts, i.e., 

• 1. Microscopic movement: we attempt to change the value of each Xi individually. We 
put x[ = Xi + dx and decide whether Xi should change or not. 

• 2. Global movement: we consider all Xj's together as a set, and consider 

^= / I- I 

aqSin ( 2Txq 

q=l 



X • = + ao + ^ aqSin \ 2Txq ^ ) , Vi = 1 ■ ■ ■ M. (4.3) 



In our simulation the total number of points in 'time' axis are chosen such that more 
points are used for small k regime. The numbers range from 7 to 16. Although this is not 
a large number, results from simulations we carried for a simple harmonic oscillator give 
pretty small errors. This gives us confidence that in the presence of random potentials the 
small number of points in the Trotter time dimension will also give us a small error in com- 
parison with the statistical error which comes from the small number of samples of random 
potentials. Since this is not a very large number of points we found it sufficient to include 
only the zeroth and first normal modes , i.e. we choose gc = 1, and then decide whether 
they move or not. The magnitude of dx and the a^'s are chosen so that the acceptance ratio 
is approximately 0.5. The parameters Oq, Oi, and the size of microscopic movement dx are 
listed in Table 1 at the end of this section. 

We have generated two kinds of random potentials: A random potential with long range 
correlation, characterized by 7 = 1/2 and a random potential with short range correla- 
tion which decay as a Gaussian. Within the framework of the variational approximation 
such a potential is equivalent to a random potential with power law correlations that are 
characterized by an index 7 = 3/2. 



For the case of long range correlation, we have generated K = 6000 uncorrelated random 
numbers hi - ■ - hx with a Gaussian distribution. We then constructed the variables Vq - ■ - Vk 
by: 

/ i K \ 

Vi = const. X l^^hj - "^hj \ (4.4) 

\i=o j=i J 

and placed them on one dimensional lattice with a lattice constant chosen to be 0.01. The 
random potential generated in this way has correlation 

(ViVj)^ = -const. X \i-j\ + C (4.5) 

where C is a constant independent of i, j. The free energy depends trivially on C but 
((x^))^ is independent of C as has been mentioned in equation (2) of ref [|1^. In this way a 
random potential with long range correlation is generated. These 6000 numbers have long 
ranged correlations with 7 = 1/2. 

We now discuss the procedure to generate a random potential with short range correla- 
tions. From the fact that random potentials with correlation 

{V{xi)V{x2))ji cc exp(-^ X (xi - xs)^) (4.6) 
in configuration space have correlations in momentum space with the following form 

{VMr (X 6ih + h) X exp(-A;iV2a), (4.7) 
we generated random numbers with a proper distribution in momentum space and fast 



Fourier transformed them by a standard Fortran subroutine p5[ to get random numbers 
with a Gaussian correlation. We put 4000 of them on lattice sites with a lattice constant 
of 0.005. These constitute a random potential with short (Gaussian) correlations. In both 
long and short-ranged disorder, we discretize the x direction such that the lattice is about 
2 orders of magnitude smaller than ((x^))^. 

Another difficulty is that in the low temperature regime the relaxation time of the system 
is very long because the phase space of the system is separated by high free energy barriers 
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and it is difficult to get reliable results within a reasonable computer time. The standard 
approach of doing a simulation in these systems is the simulated annealing approach, and 
in addition to using this method we also used a modified version of the global movement 
algorithm to speed up the dynamics. 

We started our simulation from the high temperature regime (T = 3Tc(oo) for long range 
correlation case and T = 2Tc(oo) for short range correlaion case). We then lowered the tem- 
perature in steps of ITc(oo) when T > Tc (oo). For T < Tc (oo) we lowered the temperature 
by 5T = O.lTc(oo) at a time and performed thermal averages for every 0.2Tc(oo). The lowest 
temperature of our simulation was set to be 0.2Tc(oo). In addition, we attempted 2 different 
size (zeroth mode) global movements in each sweep, of magnitudes Oq and Og respectively. 
Generally we chose one of these parameters, say oq, to be much bigger than the other (but 
such the acceptance rate will not fall below 0.01 ) in order that the particle will get a chance 
to occasionally escape from deep wells which correspond to metastable minima. We have 
performed the simulation for different given random potentials with thermalizaion sweeps 
ranging from 30,000 to 200,000 and averaged over 200,000 to 2,000,000 MC steps and only 
very small difference have been found. When taking data we have chosen the number of 
thermalizaion sweeps to be 50,000 for all cases and averaged over 200,000 to 400,000 steps 
with more steps for lower temperature. 

Table 1 gives an example for the parameters we have chosen for the simulation in the 
case of long ranged correlated potential. 



V. RESULTS AND DISCUSSION 

The mean square displacement of the particle as well as the sample to sample fluctuations 
of the mean square displacement are calculated in our simulation, and all data points are 
averaged over 2000 samples of the random potential. 

For the classical case, previous numerical results for the long range case, were reported by 
r3| for the long ranged case and also by for both the short and long ranged cases. They 
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did not perform Monte Carlo simulation, but used the fact that in the classical case instead 
of a path integral one has to evaluate a simple integral for each realization of the random 
potential. For large k{= 100), where the results should reduce to the classical case, we have 
checked our results for the mean square displacement against their's, and the agreement is 
quite good, taking into account the fact that we have averaged over 2000 realizations of the 



disordered versus 10000 in JTsland 40000 in [|l^]. In doing a PIMC we did not have enough 



computer time to average over a larger number of realizations. 

In Fig. 4 we show the results of the Monte Carlo simulations for the mean square 
displacement, for the case of long ranged correlated potential. For comparison we show the 
results obtain from the analytical solution reported in I. As in the classical case we observe 
that when quantum effects are turned on, the RS symmetric solution gives an unphysical 
result below Tc{k,), down to a certain critical k. The RSB solution gives rise to a flat 
behavior of ((x^))^ below Tc{k,). The actual results of the QMC show that the function is 
indeed monotonic, but no sharp transition is observed, and it continues to decrease at all 
temperatures. A sharp transition is only expected at = oo, whereas the simulation has 
been carried out at A^ = 1, where A^ is the number of spatial dimensions. The variational 
approximation also gives rise to a sharp transition at all dimensions, much like the large- A^ 
result. As k decreases, tunneling increases and the glassines of the system decreases as the 
particle is able to tunnel across potential barriers. This is evident in the analytical solution 
by the decrease of Tc{k) with decreasing k, until there is no longer any transition. 

In this simulation we found that the statistical errors are dominated by the sample to 
sample fluctuations of (x^) (The discussion of this quantity is in the next paragraph.). For 
this reason, the error bars in Fig. 4 are given by 

2\l/2 



±2(^(^x^)^x^))^-((x^))^j (number of samples), (5.1) 

with the number of samples being 2000 in our case. We give only the error bars for the 
case of K = 100 and k = 0.1 otherwise the figure would become too messy. The error bars 
increase with increasing k. Also, in Table 2 we give the value of {{'x'^))j^ for the first and 
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second 1000 samples for k,—0.2, 0.3, and 1.0. We found that all the values of ((x^))^ are 
within the error bars for those values of k which is very reasonable. 

What is the signature of this effect in the simulations? In order to see this effect we 
measured the sample to sample fluctuations of (x^) . This is a direct measure for the glassiness 
of the system. In Fig. 5 we depict the sample to sample fluctuation of (x^) and observe that 
for small enough k the function becomes flat with decreasing temperature which signals the 
fact that quantum effects wipe out the glassy behavior. This flgure is to be compared with 
the Fig. 2 obtained from the variational approximation. 

We also have to notice that the large sample to sample fluctuation of (x^) in the low 
temperature regime also indicates that the error bars on our graphs for (x^) are relatively 
large in the glassy phase. In fact the uncertainty of our calculation from only 2000 samples in 
not enough to make a highly precise quantitative description of the behavior of this system. 
Nevertheless, for a qualitative study of this system at low temperature our simulation gives 
a clear picture. 

Similar results were obtained for the case of short ranged correlated potential. In Fig. 
6 we display the PIMC results for the mean square displacement together with the results 
obtained in I from the analytical calculation for the RS solution. The analytical solution has 
been obtain for power correlations with 7 = 3/2, but from the variational approximation 
any faster falling correlation should give similar results. 

In Fig. 7 we show the sample to sample fluctuations of (x^) for the short ranged case, 
as obtained from the QMC simulations to be compared with Fig. 3 obtained the analytical 
variational calculation. Again we observe the reduction of the glassines of the system with 
the increase of the quantum effects. 

To conclude, we observe the relatively good agreement between the QMC simulation 
and the results of the variational calculation. Some of the deviations are of course due to 
the variational approximation, but they are also due to the fact that in the simulation we 
averaged only over 2000 realizations. The results show that the sample to sample fluctuations 
of the mean square deviation of the particle from the origin are a good measure for the 
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glossiness of the system which decreases with increasing quantum effects (increase of h'^/m). 
Thus the quahtative similarity with the phase diagram of the quantum spin glass in a 
transverse field is established. 
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APPENDIX A: 



For the temperature range where the RS solution is vahd, it has been shown in Appendix 
B of I how to obtain the numerical solution of the self consistent equations (4.15) and (4.16) 
of I. When replica symmetry is broken, in the case of long-ranged correlated potential, the 
equations for self energy matrix are given in eqs. (5.3)-(5.8) of I. The solution of equation 
(5.6) of I is obtained analytically there, and is presented in equations (5.9)-(5.11) which are 
reproduced below. 

The off diagonal elements parametrized by the Parisi variable z are given by: 



~s{z) = < 



\Az\ ^ < z < zx 



\Az^ Zi < Z < Z2 

\Az\ Z2 < z <1 



with 



A = (2/3) V/?' (Al) 
3 
T 



and Z2 is the solution of the equation 



Agrz^^ + Z2 - ^^0 . (A3) 



where 



with 



aR{P,m,iJ,,g) ^ ao + bo{P,K,iJ,,g), (A4) 



bo{f3, K,ii,g)^^Yl 2^^ ~ ( V (AS) 



We integrate s over z to get Jq dz s{z) and substitute in eq. (5.7) of I: 

ri 2 rU , / 2 1 — p-*"^'? \ 

saiu;) ^ dz s{z) -U dC{l - e^^ /' 4 ^ ' ~ ( >^ (^6) 

Jo hJo \p Jt^q mw'^ + A* - Sdii^') I 
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Then we have a set of nonhnear equations for Sd(c(j)'s which are solved for up to 20 non- 
zero Matsubara frequencies (/ = —10 • • ■ 10, in eq. |3.4| ). The solution is then used (utilizing 
eq.(5.4) of I) to calculate the sample to sample fluctuation of (x^) as given in eq. p.lO| , and 
is depicted in Fig. 2. 

For the case of short-ranged correlation, we seek a solution with one step replica sym- 
metry breaking, i.e.. 



s[z 



Sq < Z < Zc 

Si Zc < Z < 1 



and hence 



s]{z) = z s{z) - dz s{z). (A7) 







is given by 



Q < Z < Zc 
Zc < Z < 1 



where S = Zc{si — so). The breaking point and the order parameters sq, si are found 
from maximizing the variational free energy. The equations are similar to the classical case 



15| , the difference between the classical and the quantum case enters again through the 
renormalization of the constant Oq ^ = cto + which enters the correlation function of 
the random potential. 

We thus obtain a set of self-consistent equations for the Sd(u;)'s and also s{z) (or sq 
and Si) at the same time and the Newton- Raphson method can be applied. When replica 
symmetry is not broken we simply get Uc = 1. The solution is then used (utilizing eq.(5.4) of 
I) to calculate the sample to sample fluctuation of (x^) as given in eq. |3.10| , and is depicted 
in Fig. 3. 
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Table Caption: 

• Table 1: Example of different parameters in the QMC for the long ranged correlated 
potential. 

• Table 2: Example of ((x^))^ from two sets of 1000 reahzations of the long ranged 
correlated random potential in the QMC. 
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Figure Captions 



• Fig 1. Schematic phase diagram of a quantum particle in a random potential plus an 
harmonic potential, as obtained from the variational calculation derived in I. 

• Fig 2. Plot of sample to sample fluctuation of ((x^))^ for long-ranged correlated dis- 
order from numerical solution. Dashed lines are solutions assuming replica symmetry. 
Solid lines are RSB solutions. From top to bottom: k — 100, 1.0, 0.3, 0.2, 0.1. 

• Fig 3. Plot of sample to sample fluctuation of (x^) for short-ranged correlated disorder 
from the numerical solution. Dashed lines are solutions assuming replica symmetry. 
Solid hues are 1-step RSB solutions. From top to bottom: k =100, 6, 3, 2, and 1. 

• Fig 4. Plot of ((x^))^ vs. T/Tc (oo) for long-ranged correlated disorder. Solid hues are 
numerical solutions, obtained in I, assuming replica symmetry. From top to bottom: 
n =0.1, 0.2, 0.3, 1.0, 100. Data points are Monte Carlo simulation for n = 100 ( 
circles), k — 1.0 (squares), k — 0.3 (diamonds), k — 0.2 (up triangles), and k — 0.1 
(down triangles). Each point is averaged over 2000 samples. Error bars indicate the 
statistical errors for the cases of «; = 100 (dotted line) and = 0.1 (solid line). 

• Fig 5. Plot of sample to sample fluctuation of (x^) for long-ranged correlated disorder 
from Monte-Carlo simulation. Data points are for n = 100 ( circles), n — 1.0 (squares), 
K — 0.3 (diamonds), k — 0.2 (up triangles), and k — 0.1 (down triangles). Each point 
is averaged over 2000 samples. 

• Fig 6. Plot of ((x^))^ vs. T/Tc{oo) for short-ranged correlated disorder. Sohd hues are 

numerical solutions, obtained in I, assuming replica symmetry. From top to bottom: 
2, 3, 6, 100. Data points are Monte Carlo simulation for k — 100 (circles), k — Q 
(squares) , k = 3 (diamonds) , k — 2 (up triangles) , k — 1 (down triangles) . Each point 
is averaged over 2000 samples. Error bars indicate the statistical errors for the cases 
oi K, — 100 (dotted hue) and k — 1 (sohd line). 
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• Fig 7. Plot of sample to sample fluctuation of (x^) for short-ranged correlated disorder 
from Monte Carlo simulation. Data points are for k = 100 (circles), k = 6 (squares), 
K = 3 (diamonds), k = 2 (up triangles), k = 1 (down triangles). Each point is averaged 
over 2000 samples. 
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